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ON  HYPERSONIC  BLUNT  BODY  FLOW  FIELDS  OBTAINED  WITH  A 
TIME-DEPENDENT  TECHNIQUE 


Prepared  by: 

John  D.  Anderson,  Jr. 
Lorenzo  M.  Albacete 
Allen  E.  Winkelmann 


ABSTRACT:  New  results  are  presented  for  inviscid,  supersonic  and 
hypersonic  blunt-body  flow  fields  obtained  with  a  numerical  time- 
dependent  method  patterned  after  that  of  Moretti  and  Abbett.  In 
addition,  important  comments  are  made  with  regard  to  the  physical 
and  numerical  nature  of  the  method.  Specifically,  numerical  results 
are  presented  for  two-dimensional  and  axisymmetrlc  parabolic  and 
cubic  blunt  bodies  as  well  as  blunted  wedges  and  cones;  these  results 
are  presented  for  zero  degrees  angle  of  attack  and  for  a  caiorically 
perfect  gas  with  y  1.4.  The  numerical  results  are  compared  with 
other  existing  theoretical  and  experimental  data.  Also,  the  effects 
of  initial  conditions  and  boundary  conditions  are  systematically 
examined  with  regard  to  the  convergence  of  the  time-dependant  numeri¬ 
cal  solutions,  and  the  point  is  made  that  the  initial  conditions  can 
not  be  completely  arbitrary.  Finally,  in  order  to  learn  more  about 
the  performance  of  the  time-dependent  method,  a  numerical  experiment 
is  conducted  to  examine  the  unsteady  propagation  and  region  of 
influence  of  a  slight  pressure  disturbance  introduced  at  a  point  on 
the  surface  of  a  blunt  body. 


U.  S.  NAVAL  ORDNANCE  LABORATORY 
WHITE  OAK,  MARYLAND 


NOLTR  68-129 


21  August  1968 


ON  HYPERSONIC  BLUNT -BODY  FLOW  FIELDS  OBTAINED  WITH  A 
TIME-PEP END ENT  TECHNIQUE 
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perfect,  supersonic  and  hypersonic  blunt-body  flow  fields 
obtained  with  a  tine-dependent  nunurieal  method.  In  addition, 
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LI3T  OF  SYMBOLS 

speed  o f  sound 

a  4  aaa  r* -*  J  *.  M  J..  a  .  ^  > 

—  wiAjy  jiii  vne  x,  y  place 

dummy  variable 
characteristic  length 
Mach  nuaber 

pressure  (divided  by  free-stream  pressure) 
in  p 

radius  of  cylinder,  sphere,  or  paraboloid 


*  is^oftK'bSr* of  th,,  *’ y  pi*°*;  u“°  *“•-« 

t  nondlcensional  time  in  physical  coordinate  system 
T  nondimensional  time  iu  transfers-*  coordinate  system 
u  nondimensional  velocity  in  *  direction  in  a,  y  plane 

v  nondimensional  velocity  in  y  direction  in  x,  y  plane 

W  velocity  of  the  shock  in  x  direction  in  x,  y  plane 

x  abscissa  in  a  Cartesian  frai^e  in  the  physical  plane 

y  ordinate  in  a  Cartesian  frame  in  the  physical  plane 

Y  ordiniit#  In  a  Pa rtnc<«n  #«*«**<*  ■* »  >  a - „  * 

u auMunea  plane 

Cartesian  reference  frame  along  shock  or  body 
6  s  -  b 

Y  ratio  of  specific  heats 

P  density  (divided  by  free-stream  density) 

6  ancle  between  tlhirant  tn  ..J  — _ A  _ 

CTw. - -■«*■  WI*WV«  »UU  WV  QAlt f 

plane3**  *  CartesiBn  reference  frame  in  the  transformed 

♦  In  p  -  In  p 

The  subscript  -  indicates,  as  usual,  free-stream  conditions. 
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OH  HYPERSONIC  RLUNT-BQDY  FLOW  FIELDS  OBTAINED  WITH  A 
TIME-DEPENDENT  TECHNIQUE 

by 

John  D.  Anderson,  Jr. 

Lorenzo  M.  Albacete 
Allen  F.  SiakelsHum 

INTRODUCTION 

Experimental  and  analytical  Investigations  of  supersonic  and 
hypersonic  flow  fields  about  blunt  bodies  hJ*'  tcuu  prevalent  In 
aarodjmuxc  research  for  tne  past  two  decades.  This  is  in  part 
due  to  the  practical  use  of  blunt  bodies  for  manned  and  unmanned 
re-sntry  vehicles,  as  well  as  to  the  intellectually  challenging 
mathematical  nature  of  the  mixed  subsonic-supersonic  flow  field 
which  Is  characteristic  of  such  bodies.  However,  as  pointed  out 
lu  reference  (1),  it  is  remarkable  that  the  present  state-of-the- 
art  has  still  not  progressed  to  the  point  where  one  specific  ana¬ 
lytical  technique  for  solving  the  blunt-body  problem  can  be  agreed 
upon  as  being  markedly  superior  to  all  others.  On  the  other  hand, 
a  relatively  new  technique  has  appeared  in  the  recent  literature 
which  appears  to  have  the  potential  for  such  pre-eminence  among 
inviscld  blunt-body  solutions;  namely,  the  time-dependent  approach 
much  as  employed  by  Bohachevsky  et  al  (refs.  (2)  and  (3)),  and  more 
recently  by  Moretti  and  Abbett  (ret.  (4)),  and  Morettl  and  Bllech 
(rsf.  (5)).  The  present  report  supplements  the  above  investigations 
by  further  examination  and  application  of  the  time-dependent  method 
to  various  two-dimensional  and  axlsymmetrlc  blunt  bodies. 

It  Is  not  within  the  scope  of  the  present  report  to  give  a 
detailed  review  of  the  various  analytical  approaches  to  the  blunt- 
body  problem;  indeed,  excellent  descriptions  of  the  significant 
theories  are  readily  available  in  the  literature.  (For  example, 
ace  references  (1),  (6),  and  (71)  Instead,  it  is  sufficient  and 
pertinent  to  describe  briefly  the  philosophy  of  the  time-dependent 
method,  after  which  the  specific.  purpose^  of  tho  ^resent  iQpur i 
will  be  outlined. 

The  time-dependent  blunt-body  approach  entails  the  solution 
of  the  unsteady  partial  differential  equations  which  physically 
describe  the  transient  flow  field;  i.e.,  the  equations  contain 
time  derivatives  of  the  local  flow  properties.  Given  a  fixed 
body  shape  and  free-etreaa  conditions,  the  solution  begins  by 
assuming  (somswhat,  but  not  completely)  arbitrary  initial  values 
for  the  flow  field  variables  and  shock-wave  shape  and  location. 

The  solution  then  proceeds,  using  the  unsteady  equations,  to  obtain 
new  values  of  the  flow  field  variables  and  shock-wave  shape  in  steps 
of  tlms.  This  solution  is  carried  cut  numerically  by  a  finite- 
difference  scheme  using  a  high-speed  digital  computer.  Following 
a  natural  relaxation  process,  the  solution  asymptotically  approaches  s 
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steady  state  for  lar^e  values  of  time.  For  the  present  investigation, 
this  steady  state  is  the  desired  result;  the  use  of  the  time- dependent 
equations  is  just  a  means  to  an  end.  However,  the  transient  numerical 
variation  of  the  flow  field  variables  during  tneir  approach  to  the 
steady  state  should  approximate  the  actual  physical  relaxation  asso¬ 
ciated  with  a  flow  winch  is  artificially  started  with  the  assume! 
initial  conditions,  (How  well  thisr  transient  process  is  approximated 
depends  on  how  well  the  time-dependent  difference  equations  approxi¬ 
mate  me  original  governing  partial  differential  equations.)  "therefore, 
the  t iiae-dependen t  approach  is  also  useful  for  the  analyses  of 
inherently  transient  problems;  for  example,  the  establishment  of  flow 
ah-*ut  models  ir  -  tuoe ,  and  the  interaction  of  a  blast  wave 

with  a  re-entry  vehicle  (ref.  (8)).  Never theless ,  from  the  point  of 
view  of  the  steady  supersonic  blunt-body  problem,  the  mathematical 
advantage  obtained  by  use  of  the  time-dependent  method  is  that  the 
partial  differential  equations  are  hyperbolic  in  time;  this  is  in 
contrast  to  the  steady  flow  equations  which  are  elliptic  in  the  sub¬ 
sonic  region  and  hyperbolic  in  the  supersonic  region.  This  hyperbolic 
nature  of  the  time-dependent  solution  allows  the  use  of  a  finite- 
difference  approach  starting  with  assumed  initial  conditions  a,.  a 
given  time.  A.  considerably  more  detailed  description  of  the  time- 
dependent  method  can  be  found  in  references  (2)  to  (5),  and  therefore, 
no  further  elaboration  will  be  given  here.  However,  for  the  purposes 
of  establishing  nomenclature  and  for  making  this  report  somewhat 
self-contained ,  the  pertinent  equations  and  numerical  solution  employed 
in  the  present  investigation  will  be  outlined  in  a  subsequent  section. 

The  present  report  gives  new  results  for  blunt-body  flow  fields 
obtained  with  a  time-dependent  method  patterned  after  that  of  Moretti 
and  Abbett  (ref.  (4)),  and  also  makes  some  important  comments  with 
regard  to  the  physical  and  numerical  nature  of  the  method. 

Specifically,  the  main  purposes  of  the  present  report  are  fourfold; 

(1)  To  document  numerical  results  for  two-dimensional  RDd 
axlsymmetr ic ,  supertonic  and  hypersonic  flow  fields  about  various 
blunt  body  shapes;  namely  analytical  shapes  (parabolic  and  cubic)  as 
well  as  blunted  wedges  and  blunted  cones.  Results  are  presented  for 
zero  degrees  angle  of  attack  and  for  a  calorically  perfect  gas  with 
Y  -  1.4. 

(2)  To  compare  numerical  results  with  other  existing  theo¬ 
retical  and  experimental  data. 

(3)  To  systematically  describe  the  effects  of  different 
assumed  initial  conditions  and  differently  obtained  boundary  conditions 
on  the  convergence  of  the  time-dependent  numerical  solutions. 

(4)  To  describe  the  results  of  a  numerical  experiment 
conducted  to  answer  the  following  question:  If  a  slight  pressure 
disturbance  is  introduced  at  a  point  in  the  supersonic  portion  of 
the  blunt-body  flow  field,  what  is  its  region  of  influence  and  how 
clearly  is  it  defined  by  the  time-depeuden v  analysis  employed  in 
the  present  investigation?  The  answer  to  this  question  has  signifi¬ 
cance  with  regard  to  the  type  of  body  configuration  which  _an  be 
realistically  treated  by  the  above  analysis. 
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ANALYSIS 

Basic  Equations  and  Transf ormatlone 

Since  our  method  of  solution  is  patterned  after  Morettl  and 
Abbett's  annlveis,  described  in  reference  (4),  we  eive  here  only  a 
synopsis  of  the  procedure,  intended  primarily  to  permit  discussion 
of  our  results.  The  reader  is  thus  referred  to  reference  (4)  for 
the  details. 


We  are  considering  the  inviscid  flow  arcund  a  body,  as  shown 
In  figure  la.  The  body  BC  (eitner  two-dimensional  or  axisymmetric) 
is  assumed  symmetric  with  respect  to  the  x  axis;  AD  is  an  arc  of 
the  unknown  shock  wave;  the  location  of  C  is  such  that  the  flow  is 
always  supersonic  along  CD;  the  flow  at  infinity  is  uniform,  super¬ 
sonic,  and  parallel  to  the  x  axis. 


The  governing  equations  are: 


Continuity 

If 
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dy 
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In  equation  (1),  K  -  0  for  the  two-dimensional  configuration, 
and  K  -  1  for  the  axlsymmetric  case.  In  equation  (4),  *  -  In  p  - 
Yka  p.  The  energy  equation  thus  corresponds  to 


for  a  calorlcally  perfect  gas.  The  additional  symbolism  is 
summarized  in  the  list  of  symbols  at  the  beginning  of  this  report. 

C  *  ^  ,  y  ■  y;  t  ■  t  (5) 


where  <5  -  s-b,  the  local  shock  detachment  distance.  The  definitions 
(5)  transform  the  "oblique"  region  of  interest  ABCD  into  a  rectangular 
region,  as  shown  in  figure  lb. 
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In  addition,  the  following  variables  are  defined: 


C  s  (C-D  -  C  ctn  0 
P  »  In  p 
R  ■  in  p 

where  W  is  the  velocity  of  the  shod:  wave. 


B 


!~U  -  tf£  f  vC  I  /  6 
L  4 


(6) 


All  the  variables  ere  uon-dimensionailzed  as  follows:  p  and  p 
are  divided  by  their  f  recs-streans  values;  the  velocities  are  divided 
by  [p  i;  the  leap  the  are  divided  by  a  churacteristic  length  L; 

aid  thus  the  non-dimensional  time  is  the  dimensional  time  divided  by 
L  /  Sv^fT-  . 


All  of  then©  operations  result  in  the  following  non-dimensional, 
transformed  equations: 


Continuity 
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el 

-  Tb 
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A# 
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•a 

(10) 

Now  in  the  spirit  of  reference  (4),  the  flo*  field  is  divided 
into  three  parts:  "’ho  shock  itself  (AD),  the  body  (BC) ,  and  the 
’'inner  points"  insiue  the  rectangle  A3CD. 


Inner  Points  and  Finite  Differences 


We  consider  first  the  inner  points.  As  indicated  in  the 
introduction,  the  use  of  the  unsteady  equations  (hjperhoH.c  ia  time) 
enables  us  to  employ  a  finite-dil i erence  approach  j  t&rting  from  the 
initial  conditions  at  an  initial  time.  The  particular  finite- 
dif/ erence  form  is  a  matter  of  choice.  In  this  analysis,  after 
Moretti  ard  Abbett,  we  use,  for  any  flow  variable  g: 


g  (T  +  AT) 


g  ('i)  + 


AT 


(T) 


(AT)2 

2 


(ID 


That  is,  given  initial  conditions  for  each  variable,  equation  (!!) 
enables  us  to  obtain  its  value  at  a  later  time.  These  values  serve 
in  turn  as  "initial  conditions"  fer  the  next  time  step,  and  the 
process  continues  until  the  unsteady  terms  are  negligible.  This  will 
then  be  the  final  solution  for  the  steady  problem,  approached 
asymptotically  with  time. 
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All  the  terms  in  the  right-hand  side  of  (11)  are  known:  j(T) 

Is  known  from  the  prevlors  time  step  and  the  first  time  derivatives 
are  given  by  the  governing  equations  (7)  and  (10) .  In  these  equations, 
the  »pace  derivatives  are  computed  by  central  finite-differences  at 
each  pcint  in  a  mesh  which  subdivides  the  Inner  region.  The  second 
time  derivatives  are  obtained  by  differentiating  equations  (7) 
thx  ough  (10) .  Symmetry  is  assumed  at  the  centerline  for  the  calcu¬ 
lation  of  the  space  derivatives. 

For  the  axlsymuetric  case,  there  will  be  indeterminate  terms 
at  the  centerline  of  the  form  0/0.  These  may  be  calculated  by  the 
use  of  L'Hospital  rule. 

In  these  calculations,  the  step  size  AT  is  evaluated  as  the 
minimum  of  A/1.5  a  (M+l)  where  A  is  the  smallest  interval  between 
Ay  and  Ax,  a  is  the  speed  of  sound,  and  M  the  Mach  number.  This 
satisfies  the  Courant-Friedrichs-Lewy  stability  criterion. 


Finally,  note  that  in  obtaining  the  second  time  derivatives  by 

An 

differentiating  equations  (7)  to  (10) ,  the  term  —  will  be  encoun- 

dW 

tered.  From  equation  (6;  this  implies  that  —  must  be  known.  This 
term  is  computed  as  follows: 


/AT  (12) 

This  means  that  every  time  the  inner  points  are  calculated,  the  value 
of  W(T+AT),  (that  is  the  "new"  value  of  W) ^  must  be  known.  For  this 
reason,  the  shock  points  must  be  computed  first,  before  the  interior 
points,  at  the  beginning  of  each  step  in  time. 


dW 

dT 


-  W 


(T  +  AT)  -  W  (T) 


'] 


Shock,  Bouy  Points,  and  Characteristics 

The  values  of  the  flow  variables  behind  the  shock  at  T  +  AT 
can  be  obtained  from  the  Rankine-Hugoniot  relations  for  a  moving 
shock.  Bui  this  means  that  a  value  tor  W (T+AT)  must  be  assumed^ 
since  it  1b  not  known  at  the  beginning  of  each  step  in  the  compu¬ 
tation.  However,  the  values  for  the  shock  points  can  be  obtained 
in  another  way,  compared  to  those  obtained  from  the  shock  relations 
with  the  assumed  W,  and  matched  by  an  Iterative  process  until  they 
do  agree.  This  will  give  the  true  value  of  W(T+AT)  as  well  as  the 
lest  of  the  unknowns.  In  the  analysis  of  Moretti  and  Abbett,  this 
"other  way"  is  by  means  of  the  characteristic  equations  fcr  the  set 
of  governing  relations  written  for  a  (§,  r),  t)  frame,  where  §  and  rj 
are  Cartesian  coordinates  normal  and  tangential  to  the  shock  respec¬ 
tively.  The  assumption  is  made  in  obtaining  the  characteristics 
that  the  governing  equations  can  be  written  as  quasi-one-dimensional , 
modified  by  "forcing  terms"  containing  derivatives  in  the  tangential 
direction.  That  is,  the  characteristics  are  found  in  the  (s,  t) 
plane.  For  the  axi symmetric  case,  the  same  assumptions  are  made. 
However,  there  exist  In  the  characteristic  equation  ar.  added  "forcing 
term"  arising  from  the  extra  term  in  the  continuity  equation.  This 
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la  ~  where  W  Is  the  component  of  the  velocity  in  the  w  direction 

psrpcrdiculir  to  the  {?,  ")  p ir  ae .  This  term  has  a  value  of  v/y,  as 
can  be  determined  from  geometric  considerations.  At  any  rate,  the 
characteristic  equation  is  still  quasl-one-dimensional ,  that  is,  the 
characteristics  are  st^.11  found  in  the  (5,  t)  plane.  The  details  of 
the  iteration  procedure,  the  determination  of  the  characteristic 
lines,  and  the  integration  of  the  compatibility  equation,  are  given 
in.  references  (4)  and  (5). 

A  similar  computation  is  performed  at  the  body.  The  analysis 
is  simpler,  since  the  body  wall  is  fixed. 

Finally,  the  values  for  the  points  on  the  upper  boundary  are 
extrapolated  linearly  from  the  inner  values. 

All  the  points  in  the  flow  field  are  thus  account*  *'<  for. 


RESULTS 

In  general,  the  following  results  serve  to  complement  those  of 
Moretti  and  Abbett  (ref-  (4)),  and  also  lead  to  significant  comments 
with  regard  to  the  nature  and  practical  use  of  the  time-dependent 
approach  of  reference  (4) .  In  addition,  interesting  comparisons  are 
made  with  existing  theoretical  and  empirical  results,  Including  the 
useful  correlations  by  Billig  (ref.  (9))  for  shock-wave  shape  and 
detachment  distance.  In  particular,  flow-field  results  have  been 
obtained  for  the  following  blunt  bodies  at  zero  degrees  angle  of  attack: 

(1)  Parabolic  cylinder  at  Mb  ■  4  and  8,  with  the  cross  section 
given  by 

b  -  0.769  y2  -  1.0 

(2)  Cubic  cylinder  at  -  4,  with  the  cross  section  given  by 

b  -  0.427  y3  -  1.0 

(2  /  £  1  uu  t £ u  irou^c  x  .  u  Uf  ,  V* U U ra  xo  b  Xiig  XJ  JX  It  X  X XlUX  Xb-lt  A  U O &C 

of  unit  radius,  a  transition  section,  and  a  wedge  at  an  angle  of 
14  degrees  with  the  horizontal.  As  explained  later,  the  small 
transition  section  is  required  to  smooth  the  otherwise  discontinuous 
curvature  at  the  junction  point  of  the  circular  cylinder  and  the 
wedge.  This  transition  section  is  graphically  determined,  and  its 
deviation  from  the  purely  cylinder-wedge  shape  is  small. 

(4)  Paraboloid  of  revolution  at  M  -  4  and  8,  with  the  cross 
section  given  by 

b  «  0.769  r2  -  1.0 

(5)  Sphere  cone  at  *  4  and  8,  consisting  of  a  spherical 
nose  of  unit  radius,  a  transition  section  determined  graphically, 
and  a  cone  afterbody  with  a  14  degree  half-angle. 
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All  of  the  following  results  have  been  obtained  for  a  calorically 
perfect  gas  with  y  -  1.4.  In  addition,  the  8ame  grid  size,  wan 
used  for  ail  the  numerical  calculations;  namely,  a  rectangle  (see 
figure  lb)  with  C  from  0  to  1  and  Y  from  0  to  2,  with  equal  increments 
of  1/7  along  both  the  C  and  V  axis  respectively. 

The  results  will  be  given  in  four  sections.  First,  new  results 
for  several  two-dimensional  blunt-body  flow  fields  will  be  presented 
and  compared  with  other  available  theories.  Then,  significant  results 
pertaining  to  the  propagation  of  a  pressure  disturbance  will  be  shown 
and  interpreted.  Next,  the  latitude  within  which  the  initial  condi¬ 
tions  can  be  "arbitrarily"  assumed  will  be  examined  and  discussed 
and  the  affect  of  differently  obtained  boundary  conditions  will  be 
presented.  Finally,  the  results  for  the  axisymmetric  case  will  be 
presented  and  compared  with  some  existing  theoretical  and  experimental 
data, 

Two-Dimensional  Flow-Field  Results 

To  begin  with,  typical  paths  which  the  time-dependent  solution 
follows  to  a  steady  state  are  shown  in  figures  2  and  3  for  the  case 
of  the  parabolic  cylinder  at  Mm  “  4.  Figure  2  shows  the  transient 
behavior  of  W  on  the  centerline,  and  figure  3  shows  the  time  variation 
of  stagnation-point  pressure.  In  both  cases,  the  assumed  initial 
value  is  the  proper  steady-state  value.  Two  points  are  noted  from 
these  figures:  (1)  the  most  extreme  variations  occur  at  early  times 
where  the  "driving  potential"  towards  the  steady  state  is  the 
strongest,  and  (2)  the  steady  state  is  rapidly  approached  (and  for 
practical  purposes,  is  achieved)  at  large  values  of  time.  In  addition 
to  the  above  results,  figure  4  shows  the  time-dependent  shock-wave 
motion  for  the  parabolic  cylinder  at  Mb  ■  4.  Again,  the  rapid 
approach  to  a  steady  state  is  obvious.  After  about  300  time  steps, 
the  final  shock  shape  and  location  are  well  defined.  (The  steady- 
state  results  presented  in  this  report  are  those  obtained  after 
500  time  steps.)  Other  indexes  are  also  used  to  measure  the  degree 
to  which  the  steady  state  is  achieved,  such  as  the  variation  of 
centerline  entropy  and  position  of  the  sonic  line.  All  of  these 
indexes  indicate  that  the  steady  state  was  achieved  (practically 
speaking)  for  the  cases  given  in  the  present  report. 

Final  steady-state  surface-pressure  distributions,  normalized 
with  respect  to  stagnation  point  pressure,  are  given  in  figure  5 
for  the  parabolic  cylinder  tt  M#  ■  4  and  8.  These  results  are 
compared  with  the  modified  Newtonian  Formula  (ref.  (10)).  Figure  5 
contains  two  anticipated  results;  namely,  that  the  modified 
Newtonian  distribution  underestimates  the  actual  pressure  distribution, 
and  that  the  actual  pressure  distribution  is  closer  to  Newtonian 
at  Ma  ■  8  than  at  M  ■  4.  This  type  of  comparison  with  modified 
Newtonian  is  typical  of  two-dimensional  blunt  bodies;  in  practice, 
pressure  distributions  over  axisymmetric  bodies  agree  more  closely 
with  modified  Newtonian  than  do  two-dimensional  pressure  distributions, 
as  will  be  shown  later. 
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The  1  ioal  «L;r*tik  sbipAs  and  sonic  linnn  >r«  «hn»n  in  f i.^ure  6 
for  the  parabolic  cylinds?  at  M  -  4  and  8.  The  decreased  shock-" 
detachment  distance  led  oownward  shift  of  the  sonic  line  are  both 
characteristic  of  an  i<r. creased  Mach  number.  In  addition,  the 
resulting  shape  of  the  uouic  line  is  proper  f  or  tw;i-ril  msnaj  r.n«  1 
blunt  bodies  at  .  (For  example,  see  reference  (1).) 

Figure  7  shows  the  steady-state  streamline  direction  field 
for  the  parabolic  cylinder  at  -  4. 

Figure  8,  which  shows  the  time  variation  of  stagnation  point 
pressure  for  Mw  -  8,  demonstrates  a  particularly  interesting  point 
when  compared  with  figun  3,  which  Is  for  M  -  4.  Starting  with 
similar  initially  assctcd  £lo’«  fields  and  sfiock  shape,  figures  3 
and  8  demonstrate  that  th®  approach  to  a  steady  state  is  approximately 
a  factor  of  4  faster  at  B()  **  8  than  at  M  -  4.  In  fact,  for  M  -  8, 
InfpQ/p^)  has  converged  to  within  0.1  percent  of  its  steady-state 
value  by  t  -  1.0;  this  i*  in  contrast  to  t  -  4.0  for  similar 
convergence  at  M  ■>  4.  Taeue  results  prompt  the  following  explanation. 
The  factor  of  4  difference  In  relaxation  time  exists  because  the 
local  sound  speed  and  ihe  local  flow  velocity  within  the  shock  layer 
are  both  increased  by  approximately  a  factor  of  two  in  going  from 
•  4  to  Mff  ■  8.  Consequently,  relative  to  the  fixed  body,  the 
velocity  of  a  weak-pi  assure  wave  propagating  downstream  is  increased 
by  a  factor  of  4,  resulting  in  faster  relaxation  (by  a  factor  of  4) 
of  the  flow  field  to  the  steady  state.  The  above  phenomena  identify 
the  following  simple  gasdynamic  mechanism.  The  rapidly  changing 
unsteady  bow  shock  wave  causes  transient  pressure  changes  immediately 
behind  the  wave.  In  turn,  the  subsequent  downstream  propagation  and 
interaction  of  these  unsteady  pressure  waves  constitute  the  primary 
physical  mechanism  by  which  the  unsteady  flow  field  relaxes  to  the 
steady  state. 

Results  obtained  for  the  cubic  cylinder  at  ■  4  are  shown 
in  figures  9  and  10,  which  illustrate  the  shock-wave  shape  and 

aurfiiOe  pi«bSUr«  ulbti  ibiitiOut  Ap*&l.u,  tiiv  biiOCk**W£.vv  is 

to  rapidly  approach  a  steady  cv&te,  and  as  before,  the  surface 
pressure  distribution  is  compared  with  modified  Newtonian. 

A  two-dircensional  configuration  of  more  practical  importance 
is  that  of  a  blunted  wedge;  there fire,  such  results  are  presented  in 
figures  11-14  for  Mw  -  8.,  As  an  .index  of  the  convergence  behavior, 
figure  11  shows  the  time  variation  oi  stagnation  point  pressure. 

Again,  the  rapid  convergence  n  th:  steady  state  is  obvious; 
has  converged  to  within  0.1  percent  of  its  steady-state  value  by 
t  “  1.0.  And  again,  the  last  relaxation  associated  with  the  Mach 
number  is  noted.  In  fact,  the  approach  to  the  steady  state  and  the 
relaxation  time  shown  in  figure  11  are  comparable  to  those  shown  in 
figure  8  for  the  parabolic  cylinder,  thus  Indicating  that  relaxation 
time  is  relatively  independent  of  fcsdy  shape. 
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^Vi*;  j  -  '^HrW 


Figure  12  presents  the  time-dependent  bow-shock  wave  motion,  as 
well  as  the  final  shock  wave  and  sonic  line.  Rapid  convergence  is 
also  apparent  in  this  figure.  In  addition,  results  obtained  from  the 
correlation  of  Billig  (ref.  (9))  for  final  shock-wave  shape  and  loca¬ 
tion  are  also  presented  in  figure  12.  For  a  cylinder-wedge,  the 
correlation  is. 

f  - 1  + 1  -  [« +  1/2  - x]  Rc c,n2<J 


where  ^ 

R 

end  ^c 
R 


0.386  exp 
1.386  exp 


(4.67/M^2) 

^l.a/CM,,-])0'75] 


(13) 


6  is  the  shock  detachment  distance.  R  is  the  cylinder  radius,  Rc  is 
the  shock  vertex  radius  of  curvature,  and  6  is  the  shock  angle 
approached  asymptotically  for  large  y;  in  the  present  case,  0  is  the 
attached  oblique  shock  angle  for  the  wedge  afterbody.  The  results 
presented  in  figure  12  show  very  good  agreement  between  Billig ’s 
correlation  and  present  exact  numerical  solution. 


The  steady-state  surface  pressure  distribution  over  the  blunted 
wedge  as  a  function  of  the  vertical  coordinate,  Y,  is  shown  In 
figure  13.  (Also  shown  is  the  extent  of  the  transition  section 
mentioned  earlier.)  Two  points  are  noted  from  this  figure:  (1)  the 
proper  supersonic  weege  pressure  is  achieved  asymptotically  and 
monotonically  far  downstream  of  the  blunt  nose;  and  (2)  the  pressure 
distribution  over  the  cylindrical  portion,  compares  favorably  with 
infinite  Mach  number  circular  cylinder  results  obtained  by  Fuller 
using  an  inverse  method  (refs.  (1)  and  (11)).  At  hypersonic  speeds, 
the  pressure  distribution  is  somewhat  insensitive  to  Mn.  Consequently, 
an  infinite  Mach  number  comparison  is  reasonably  valid? 


Figure  14  illustrates  the  sensitive  transient  behavior  of  the 
sonic  line  during  convergence  to  the  steady  state.  Here  again  two 
points  ara  of  interest:  (1)  rapid  convergence  to  a  steady-state  sonic 
line  is  oi’-vious,  and  (2)  the  movement  oi  the  sonic  point  on  the  bow 
shock  wave  is  considerably  greater  than  on  the  body;  in  fact,  the  sonic 
point  movement  or  the  body  is  confined  to  a  relati fely  narrow  region. 

Figure  15  illustrates  the  steady-state  velocity  distribution 
along  the  stagnation  streamline.  The  results  are  shown  for  the 
parabolic  cylinder  at  M  ■  4  and  8,  the  cylinder-wedge  at  M  “8, 
and  the  cubic  cylinder  at  M  -  8.  The  radii  of  curvature  a?  the 
stagnation  point  for  these  bodies  are  0,35,  1  and  Infinite, 

respectively.  In  figure  15,  U_  is  the  velocity  immediately 
behind  the  normal  portion  of  th§  shock,  and  6  is  the  shock  detach-' 
ment  distance.  The  fact  that  the  velocity  distribution  is  approxi¬ 
mately  linear  has  been  noted  in  previous  analyses  (ref.  (1)). 

However,  the  exact  nonlinear  velocity  distributions  are  identical 
(within  graphical  accuracy)  for  the  given  two-dimensional  cases  (but 
this  was  not  so  for  the  axisymmetric  configurations,  as  will  be  seen 
later)  .  Consequently,  it  appears  that  for  the  two-dimensional  cases 
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the  exact  velocity  distribution  along  the  centerline  is  relatively 
independent  of  both  Mach  number  and  body  shape.  Such  independence 
has  been  predicted  by  the  closed-form  approximate  analysis  of  Li  and 
Geiger  (ref.  (12)).  However,  the  approximate  velocity  distribution 
obtained  from  Li  and  Geiger's  formula  differs  somewhat  from  the 
present  numerical  results  ar  seen  from  figure  15. 

Ob  the  Propagation  of  a  Weak  Pressure  Disturbance 

The  following  results  are  given  in  ans**»r  to  the  question  ashed 
earlier  in  the  introduction;  namely,  if  a  slight  pressure  disturbance 
is  introduced  at  a  point  in  the  supersonic  portion  of  the  blunt-body 
flow  field,  what  is  its  region  of  i  if luence  and  how  clearly  is  it 
defined  by  the  time-dependent  analysis  employed  In  tne  present 
investigation?  Towards  this  end,  a  numerical  experiment  was  conducted 
in  which,  first,  the  steady  state  was  reasonably  achieved  (after  400 
time  steps)  for  the  flow  over  the  parabolic  cylinder  at  Ma  -  4,  and 
then  a  small  protuberance  was  extended  into  the  flow  from  the  body 
surface  at  a  point  far  downstream  of  the  sonic  point,  as  shown  in 
figure  16.  The  pressure  disturbance  caused  by  this  protuberance  was 
subsequently  observed  as  a  function  of  time.  Figures  17  and  18  show 
the  time  variation  of  the  body  surface  pressure  at  the  various  points 
identified  in  figure  16.  Note  that  the  scale  is  the  same  for  all 
curves  shown  in  figures  17  and  18.  Also  note  that  the  expanded  scale 
gives  the  exaggerated  illusion  of  a  time  varying  pressure  before  the 
time  of  the  disturbance;  in  reality,  the  steady  state  has  been 
achieved  for  all  practical  purposes  by  this  time.  These  results  show 
that  the  protuberance  indeed  causes  a  noticeable  pressure  disturbance; 
this  disturbance  was  strongest  at  the  protuberance,  and  decays  both 
downstream  and  upstream  of  the  source.  Note  also  the  growing  time 
lag  of  the  disturbance  felt  at  points  further  away  from  the  protuber¬ 
ance.  Of  most  consequence,  however,  is  the  fact  that  the  effect  of 
the  protuberance,  which  is  far  inside  the  supersonic  region  of  the 
flow  field,  is  fed  upstream,  even  to  the  extent  of  the  subsonic 
portion  of  the  flow  field  and  all  the  way  to  the  stagnation  streamline. 
This  latter  effect  is  graphically  shown  in  figure  19,  which  compares 
the  time  variations  of  centerline  wave  velocity  for  the  cases  with 
and  without  the  disturbance.  As  seen  in  this  figure,  the  disturbance 
is  obviously  felt  all  the  way  to  the  centerline. 

The  upstream  influence  described  above,  which  is  contrary  to  the 
normal  behavior  in  a  steady  supersonic  flow,  can  be  attributed  to  the 
fact  that  the  insertion  of  the  protuberance  initiates  a  new  transient 
flow  field  (ref.  (13));  consequently,  all  regions  of  the  flow  field 
can  eventually  be  influenced  by  unsteady  wave  motion.  In  addition, 
even  though  the  disturbed  pressures  shown  in  figures  17  and  18  have 
not  settled  to  a  steady  state  at  the  largest  time  considered,  the 
pressures  at  the  upstream  locations  appear  to  experience  a  "permanent 
set,"  i.e.,  the  pressures  seem  to  be  approaching  steady-state  values 
different  from  the  original  steady-state  values  before  the  distur¬ 
bance.  This  "pseudo-elliptic"  behavior  can  be  attributed  to  a 
numerical  effect;  namely,  the  internal  wave  caused  by  the  disturbance 
is  smeared  over  several  mesh  points,  and  because  a  course  grid  is 
is  employed,  the  wave  can  be  felt  some  distance  upstream.  If  at  least 
one  of  these  mesh  points  is  upstream  of  the  sonic  point  (as  in  the 
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present  case)  then  the  entire  steady-s~ate  subsonic  region  will  be 
affected.  The  above  results  have  significant  practical  implication; 
namely,  that  tbe  present  analysis  should  be  applied  to  smooth  con¬ 
figurations. 

It  must  be  emphasized  that  the  present  results  do  not  degrade 
or  detract  from  the  excellent  and  useful  method  presented  in 
reference  (4);  quite  the  contrary,  the  present  authors  have  already 
stated  their  opinion  in  the  Introduction  that  the  time-dependent 
method  (such  as  formulated  by  lioretti  and  Abbett)  may  become  pre¬ 
eminent  among  practical  blunt-body  solutions.  Instead,  tbe  present 
results  provide  some  indication  that  the  method  (in  its  present 
formulation)  should  be  restricted  to  smooth  bodies.  Finally,  it  is 
interesting  to  note  that  the  method  of  reference  (4)  does  indeed  sense 
the  presence  of  a  disturbance  or  internal  wave,  and  such  a 
be  traced  rather  coar-cly  by  a-a-uo  o£  the  peaks  in  the  pressure  and 
density  distributions. 

Experiments  on  the  Sensitivity  of  the  Blunt-Body  Solution  to  Initial 
Conditions,  Body  Shape,  and  Boundary  Conditions 

In  addition  to  the  above  numerical  experiment,  other  experiments 
were  performed  to  investigate  the  effect  of  initial  conditions,  body 
shape,  and  boundary  conditions  on  the  time-dependent  blunt-body 
solutions.  These  experiments  indicate  that  such  factors  have  a 
definite  influence  on  tns  "go"  or  "no-go"  nature  of  the  blunt-body 
computer  program.  The  assumed  initial  flow  field  (at  time  equal  to 
z-jio)  is  obtained  from  the  following:  (1)  modified  Newtonian  pressure 
distribution  along  the  body  surface,  (2)  isentropic  expansion  along 
the  body  streamline,  (3)  assumed  shock  shape  snri  detachment  distance, 
and  (4)  linear  interpolation  between  the  shock  and  the  body. 

a.  Effect  of  Initial  Pressure  Distribution 

Solutions  for  a  two-dimensional  parabolic  body  were  obtained 
using  initial  pressure  distributions  which  were  factors  of  six  above 
and  below  the  standard  solution.  In  each  case,  the  solution  converged 
to  the  same  final  results.  From  this  it  is  concluded  that  the  program 
is  not  very  sensitive  to  initial  pressure  distributions,  which  can  be 
prescribed  within  an  order  of  magnitude  of  the  final  pressure  distri¬ 
bution. 

b.  Effect  of  Shock  Shape 

Satisfactory  results  were  obtained  for  a  two-dimensional 
parabolic  body  using  both  parabolic  and  cubic  initial  shock  shapes, 
with  6  -  .37  at  M  ■  4.  However,  the  use  of  a  hyperbolic  shock  shape, 
5  =  .37,  did  not  Tead  to  a  converged  solution.  Hence,  the  initial 
shock  shape  may  i;?ve  some  influence  on  the  program,  and  it  appears 
reasonable  to  run  similar  analytic  shapes  for  both  the  body  and  the 
shock.  For  blunted  wedges  and  cones  initial  parabolic  shock  shapes 
were  satisfactory. 
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c.  Effect  of  Shock-Detachment  Distance 


The  Initial  shock-detachment  distance,  ^initial  ls  an 
Important  quantity,  and  can  be  prescribed  only  within  a  certain 
latitude.  For  example,  using  a  two-dimensional  parabolic  body  and 
initial  shock  shape  at  “  4  the  program  worked  successfully  for 
6  initial  -  0.37  and  0.49,  but  it  did  not  run  for  6  initial  -  0.55 
or  0.20.  Consequently ,  the  initial  shock  cannot  be  assumed  too  close 
or  too  far  away  from  the  body.  If  indeed  6  initial  is  too  large  or 
too  snail,  the  shock  wave  tends  to  accelerate  considerably  faster 
than  shown  on  figure  2  at  early  values  of  time  (accelerations  an 
order  of  magnitude  faster  can  occur).  In  turn,  strong  gradients  of 
the  flow-field  variables  are  produced  behind  the  highly  accelerating 
wave,  and  consequently  the  finite-difference  scheme  using  a  fixed 
course  grid  becomes  inaccurate.  This  Inaccuracy  ultimately  causes 
•ud*  —  *>!> sd  wjl  vLo  prc-r^m  to  collapse. 


d.  Effect  of  Body  Shape 

Successful  solutions  were  obtained  for  bodies  witn  constant 
and  linear  variations  of  d2b/dy2,  namely  parabolic  and  cubic  shapes, 
respectively.  However,  the  program  did  not  run  for  those  parabolic 
and  cubic  shapes  which  were  characterized  by  high  values  of  d2b/dy2 
(above  approximately  6).  In  addition,  the  program  did  not  run  for 
bodies  with  local  discontinuities  in  d2b/dy2,  namely,  purely 
hsnicylinder-wedge  or  sphere-cone  shape^.  The  program  became  unstable 
In  the  region  of  the  discontinuity  in  d“b/dy  ,  i.e.,  where  the  body 
curvature  is  discontinuous.  When  a  transition  section  was  added  to 
the  discontinuity,  succsssful  solutions  were  obtained. 


e.  Effect  of  Extrapolation  to  the  Upper  Boundary 


As  shown  in  figure  la  and  as  discussed  earlier,  the  upper 
boundary  of  the  flow  field  is  taken  as  a  horizontal  line  within 
the  supersonic  region.  At  any  given  time  step,  flow-field  variables 
along  this  boundary  are  obtained  by  linear  extrapolation  from  points 
below.  Consequently,  the  upper  boundary  values  do  not  exactly 
satisfy  the  continuity,  momentum  and  energy  equations,  and  therefore 
final  results  for  certain  portions  of  the  flow  field  may  be  com- 
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effect,  results  obtained  by  using  several  different  extrapolation 
procedures  were  compared.  These  extrapolation  procedures  are: 


(1)  linear  extrapolation 

(2)  reduced  linear  extrapolation  (extrapolated  increment 
reduced  26  percent  below  that  for  strictly  linear  extrapolation) 

(3)  Shank's  approximation  (ref.  (14)) 


We  consider  first  the  two-dimensional  case.  Preliminary  data  were 
obtained  for  a  cylinder-wedge  in  helium  (y  -  1.67)  at  -  16.  These 
data  show  that  procedures  (1)  and  (2)  give  different  surface  pressure 
results  downstream  of  the  right  running  characteristic  from  point  D, 
as  shown  schematically  in  the  diagram  below. 
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Also,  procedure  (3)  resulted  In  a  collapse  of  the  computer  program, 
hence  no  results  were  obtained.  Obviously,  for  this  case,  the  pro¬ 
cedure  for  extrapolation  to  the  upper  boundary  influences  that  part 
of  the  flow  field  downstream  of  the  right  running  characteristic 
from  D  (shaded  region  in  the  above  sketch) ,  and  this  fact  should  be 
kept  in  mind  when  interpreting  the  present  numerical  results. 

However,  this  anomalous  behavior  seems  to  be  limited  to  the  flow 
field  around  composite  shapes  (such  as  a  cylinder-wedge  and  sphere- 
cone)  ;  the  flow  field  around  purely  analytical  shapes  appears  to  be 
rather  insensitive  to  the  extrapolation  procedure,  as  will  be 
discussed  further  in  a  subsequent  section  dealing  with  axisymmetric 
results.  A  complete  resolution  and  understanding  of  this  effect  is 
deferred  for  future  investigation. 

As  an  interim  summary,  the  convergent  or  divergent  behavior 
of  the  time-dependent  blunt- body  solution  appears  to  be  somewhat 
dependent  upon  the  initial  conditions  and  body  shape,  outside  of  a 
certain  latitude.  Withiu  this  latitude,  successful  solutions  were 
obtained.  However,  it  should  be  remembered  that  a  fixed  grid  size 
(7x14)  was  employed  for  all  the  present  results.  A  finer  grid  might 
reasonably  lead  to  an  even  wider  latitude  for  successful  initial 
conditions  and  body  shapes.  Also,  it  is  apparent  that  a  substantial 
portion  of  the  supersonic  flow  field  is  affected  by  the  extrapolation 
procedure  to  the  upper  boundary. 

Axisymmetric  Flow  Field  Results 

As  for  the  two-dimensional  case,  typical  paths  which  the  three- 
dimensional  time-dependent  solution  follows  to  a  steady  state  are 
shown  in  figures  20,  21,  22,  and  23.  Figures  20  and  21  show  the  time 
variation  of  wave  velocity  at  the  centerline  for  the  paraboloid  of 
revolution  at  M  “4  and  s.  respectively.  Qualitatively,  both  figures 
are  similar  to  the  two-dimensional  case,  in  which  the  most  extreme 
variations  occur  at  early  times,  but  steady  state  is  achieved  earlier 
in  the  axisymmetric  case  due  to  the  three-dimensional  relieving  effect. 
Figures  22  and  23  show  the  time  variation  of  the  stagnation  point 
pressure  for  the  paraboloid  at  Hm  ■>  4  and  8,  respectively. 
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Figures  24  and  25  present  tbe  time-dependent  shock-wave  motion 
for  the  paraboloid  it  *  4  and  a,  respectively.  After  about  200 
time  steps,  the  final  shock  shape  and  location  are  well  defined. 

In  figures  26  and  27  we  have  plotted  the  final,  steady-state 
surface  pressure  distributions  (divided  by  the  stagnation  point 
pressure)  for  the  paraboloid  at  M  -  4  and  6,  respectively.  We  can 
see  that  the  Newtonian  pressure  distribution  is  in  better  agreement 
with  tbe  axisymmetric  case  than  with  the  two-dimensional  case  (cf. 
figure  5).  In  addition,  the  pressure  distribution  is  compared  with 
results  obtained  with  tbe  inverse  method  of  Lomax  and  Inouye  (ref.  (15)). 
(All  data  from  this  source  were  kindly  provided  by  Mr.  Robert  Thompson 
of  tbe  Naval  Ship  Research  and  Development  Center,  Washington,  D.  C.) 

Figures  :8  and  ?9  present  results  for  the  sphere-cone  case 
at  M  ■  4  and  8,  respectively.  We  can  see  again  the  close  agreement 
between  tbe  modified  Newtonian  distribution  and  the  t  re-dependent 
results.  These  figures  also  show  tbe  overexpansion  and  subsequent 
approach  to  the  cone-pressure  value,  as  indicated  by  PC/PQ,  Note 
that  tbe  abflicisaa  in  these  figures  is  ‘'s'\  tbe  distance  along  the 
surface  of  tbe  body. 

In  addition,  results  from  reference  (16)  (which  uses  a  program 
similar  to  Lonax  and  Inouye 's,  coupled  with  a  method  of  characteristics 
fiowntstream)  are  included  in  figure  29.  We  see  that  agreement  is 
good  up  to  9  —0.8,  a  "ter  which  our  pressure  distribution  falls  below 
both  tbe  Newtonian  distribution  and  that  of  reference  (16) .  We 
believe  this  to  be  due  to  tbe  presence  of  our  transition  section, 
which  begins  near  this  point.  Tbe  exact  location  of  tbe  transition 
section,  and  further  comparisons  with  other  methods,  appear  in  the 
next  two  figures.  Figures  30  and  31  present  the  same  results  as 
28  and  29  but  covering  a  smaller  region  of  the  s  coordinate  shown  in  an 
enlarged  scale.  This  was  done  to  facilitate  comparison  of  our  results 
with  existing  theoretical  and  experimental  data  for  spheres.  In 
particular,  figure  30  compares  the  present  time-dependent  calculations 
for  M  “4  with  the  inverse  methnr!  of  reference  (15) ;  the  tise- 
depenclent  results  of  Moretti  and  Abbett  (ref.  (4)),  the  analytical 
results  of  Van  Tuyl  obtained  with  rational  approximations  (ref.  (17)), 
and  the  experimental  data  of  Stallings  (ref.  (18)).  Good  agreement 
is  obtained  over  most  of  the  spherical  portion;  however,  some  effect 
in  the  region  of  the  transition  section  is  again  noticeable. 

Figure  31  shows  similar  results  for  -  6. 

Figure  32  shows  the  final  shoex  shapes  for  the  sphere-cone 
at  M  -  4  and  8.  As  in  the  case  of  the  cylinder-wedge,  excellent 
agreement  with  the  correluLiuii  ui  dIxh^  xOi-  i»he  snock  shape  is 
obtained  for  the  sphere-cone.  In  addition,  the  shock  shapes  are 
compared  to  the  ones  obtained  with  the  method  of  reference  (15)  in 
figures  33  and  34  for  Mw  -  4  and  8,  respectively.  Once  more, 
excellent  agreement  is  noted. 
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A  comparison  of  shock-detachment  distance  is  mane  with  tns 
inverse  method  oi  Lomax  and  Inouye  (ref.  (15))  in  the  iolJ owing  table. 
H  is  the  local  radius  of  curvature  at  the  stagnation  point.  The 


agreement  is 

quite  good. 

X  /tv 

V/  JTi 

Shape 

K 

Ref.  (15)) 

Present  re 

Sphere 

4 

.1753 

.1760 

8 
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Paraboloid  ,, 
b  -  0.769  y* 

4 

.1264 

.1284 

-  1.0 

8 

.098 

.100 

Paraboloid 
b  -  0.450  y 

8 

.1689 

.1687 

-  1.0 

Tne  velocity  distributions  for  the  sphere-cone  at  Um  -  4  and  8 
are  given  in  figures  35  and  36,  respectively.  The  results  are  com¬ 
pared  to  Van  Tuyl's  distributions  (ref.  (17))  and  indeed  excellent 
agreement  is  observed. 

For  the  axisymmetric  case,  figures  37  and  38  (which  correspond 
to  figure  15  for  the  two-dimensional  case)  give  the  steady-state 
velocity  distribution  along  the  stagnation  streamline.  Indeed,  the 
velocity  distribution  is  approximately  linear  as  expected, but,  unlike 
the  two-dimensional  results,  the  distributions  are  not  exactly 
independent  of  Mach  number  and  body  shape  (although  they  are  close) . 
Once  more,  the  approximate  velocity  distribution  obtained  from  Li 
and  Geiger's  formula  (ref.  (12))  differs  from  our  numerical  results. 

Figure  39  compares  our  sonic  line  for  the  -  4  paraboloid 
to  the  sonic  line  obtained  with  the  method  of  reference  (15) .  The 
agreement  is  excellent.  The  agreement  was  not  so  good,  however,  for 
the  sphere-cone  case,  probably  because  of  the  presence  of  the  transi¬ 
tion  section,  located  .indeed  in  tne  neighborhood  of  the  sonic  point. 

Finally,  figure  40  illustrates  the  effect  on  the  surface 
pressure  distribution  about  a  paraboliod  at  -  4  of  various  extra¬ 
polation  procedures  to  the  upper  boundary.  We  can  see  that  these  have 
little  effect.  On  the  other  hand,  the  sphere-cone  case  is  sensitive 
to  the  kind  of  extrapolation  procedure  used;  for  example,  a  converging 
solution  using  Shank's  approximation  could  not  be  obtained.  Taken 
in  conjunction  with  similar  results  for  the  two-dimensional  case 
(cf.  (e) ,  p.  12),  we  are  led  to  believe  that  the  flow  field  around 
composite  shapes  (i.e.,  where  there  exist  discontinuities  in  the 
curvature)  are  particularly  sensitive  to  instabilities.  But  the  flow 
field  around  purely  analytical  shapes  does  not  exhibit  this  sensi¬ 
tivity. 
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New  results  tor  calorically  perfect,  supersonic  and  hypersonic 
flow  fields  about  blunt  bodies  at  zero  degrees  angle  of  attach  have 
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and  which  prompt  some  significant  comments  with  regard  to  the  ohysical 
and  numerical  nature  of  the  present  time-dependent  method.  These 
results  lead  to  tbe  following  conclusions: 


(1)  Pressure  waves  which  originate  behind  the  initially 
unsteady,  moving  bow  snock  and  which  subsequently  propagate  and 
interact  downstream  appear  to  the  principal  physical  mechanism 
by  which  the  unsteady  flow  field  converges  to  tbe  steady  state. 
Consequently,  higher  free-stream  Mach  numbers  (and  thus  higher  values 
of  the  local  speed  of  sound  and  flow  velocities  within  the  flow 
field)  result  In  proportionally  shorter  physical  convergence  times. 


(2)  k  pressure  disturbance  caused  by  a  protuberance  which 

is  suddenly  introduced  into  the  steady  supersonic  region  can  propagate 
upstream  as  well  as  downstream,  and  can  eventually  affect  all  regions 
of  the  flow  field.  This  behavior  is  a  result  of  the  unsteady  wave 
motion  F.nd  new  transient  flow  field  initiated  by  the  disturbance. 
However,  the  disturbed  upstream  surface  pressures  appear  to  approach 
new  steady  -state  values  which  are  displaced  from  the  steady-state 
values  before  the  disturbance.  This  apparent  "pseudo-elliptic" 
behavior  is  attributed  to  the  smearing  out  of  the  wave  over  the 
rather  coarse  grid  used  in  the  finite-difference  scheme, 

(3)  Billig's  empirical  correlation  for  bow  shock  shape  and 
location  agrees  very  well  with  the  present  numerical  solutions  for 
cylinder-wedges  and  sphere-conos, 

(4)  Whether  the  numerical  solution  converges  or  diverges 
depends  to  some  extant  on  the  assumed  initial  conditions.  These 
initial  conditions,  such  as  bow  shock-wave  shape  and  location,  can 
be  arbitrarily  assumed  within  a  relatively  wide  latitude.  However, 
beyond  this  latitude,  the  initial  conditions  do  not  1* ad  to  conver¬ 
gence  Therefore,  lor  practical  numerical  reasons,  tne  initial 
conditions  can  not  be  completely  arbitrary. 

(5)  Convergent  solutions  could  not  be  obtained  about  body 
shapes  with  ft  local  discontinuity  in  the  surface  curvature.  There¬ 
fore,  the  present  numerical  program  is  restricted  to  relatively 
smooth  bodies. 

(6)  The  flow  field  around  composite  shapes  (such  as  a 
cylinder-wedge  or  spheie-coii)  is  particularly  sensitive  to  instabil¬ 
ities,  whereas  the  flow  field  around  purely  analytical  shape u  does 
not  exhibit  this  sensitivity. 

Sven  though  some  of  the  above  conclusions  delineate  restric¬ 
tions  on  the  time-dependent  method  as  employed  in  the  present  investi¬ 
gation,  emphasis  must  jc  made  that  all  existing  blunt-body  solutions 
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are  restricted  in  some  roanect 

the  time-dependent  method  has  the  advan^^o?  Tire^S”' 

-f  and  free-stream  properties  are  given)  as  well  as  being 
able  to  calculate  large  portions  of  the  supersonic  flow  field 

downstream  of  the  sonic  line  /thn  ayi  ui _ j.  • _ . 

methods  usually  encounter  diff 

invri™1*  liUGi  ’nJ"  addition,  the  time-dependent  method  does  not 
nnH^in'  tsf  ^“P1  ^cation  or  reduction  of  the  equations  of  change, 
and  in  this  respect  it  is  an  exact  solution.  Finally,  it  is 
anticipated  that  the  time-dependent  method  will  be  able  to  treat 
more  complex  body  shapes  than  the  existing  steady-state  methods. 

.  r  these  reasons,  and  from  experience  gained  during  the  present 
investigation,  the  time-dependent  solution  appears  to  have  the 
potential  for  pre-eminent  application  to  the  supersonic  and  hyper¬ 
sonic  blunt- body  problem.  yv 

Final. ly ,  since  the  completion  of  the  present  studies, 
additionai  information  concerning  the  time-dependent  blunt-body 
method  has  very  recently  been  presented  by  Moretti  (ref.  (19)) 

In  general,  reference  19  confirms  many  of  the  findings  of  the  present 
studies,  and  in  addition,  il  presents  a  large  number  of  blunt-body 
flow-field  calculations  in  a  systematic  form  for  a  range  of  free- 
stream  Mach  numbers  and  body  shapes* 
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FIG.  4  TIME  DEPENDENT  SHOCK  WAVE  MOTION,  PARABOLIC  CYLINDER,  M 


FIG.  5  SURFACE  PRE5SURE  DISTRIBUTIONS;  PARABOUC  CYLINDER 
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FIG.  8  TIME  VARIATION  OF  STAGNATION  POINT  PRESSURE, PARABOLIC  CYLINDER,  M, 


FIG.  9  TIME  DEPENDENT  SHOCK  WAVE  MOTION;  CUBIC  CYLINDER, M,^  =4 
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FIG.  10  SURFACE  PRESSURE  DISTRIBUTION;CUBtC  BODY,  =  4 


FIG.  11  TIME  VARIATION  OF  STAGNATION  POINT  PRESSURE;  BLUNTED 
WEDGc,  M  =8 
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FIG.  12  TIME  DEPENDENT  SHOCK  WAYE  MOl  ION;  BLUNTED  WEDGE, IW.  =B 
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FIG.  13  SURFACE  PRESSURE  DISTRIBUTION;  BLUNTED  WEDGE, “8.0 
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FIG.  15  VELOCITY  VARIATION  ALONG  THE 
STAGNATION  STREAMLINE 
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FIG.  16  CONFIGURATION  FOR  PROTUBERANCE  EXPERIMENT 
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FIG.  29  SURFACE  PRESSURE  DISTRIBUTION;  SPHERE-CONE,  M 
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11.  SUPPLEMENTARY  NOTES 


j  12-  SPONSORING  MIL!  tAW  Y  ACTIVITY 


New  results  are  presented  for  inviscid,  supersonic  and 
hypersonic  fclunt-body  flow  fields  obtained  with  a  numerical  time- 
dependent  method  patterned  after  that  of  Moretti  and  Abbott.  In 
addition,  important  comments  ave  made  with  regard  to  the  physical 
and  numerical  nature  of  the  method.  Specifically,  numerical  results 
are  presented  for  two-dimensional  and  axisvmroetrie  parabolic  and 
cubic  blunt  bodies  as  well  as  blunted  wedges  and  cones;  these  results 
are  presented  for  zero  degrees  angle  of  attack  and  for  a  caloric  lly 
perfect  gas  with  y  -  1.4.  The  numerical  results  are  compared  with 
other  existing  theoretical  and  experimental  data.  Also,  the  effects 
of  initial  conditions  and  boundary  conditions  are  systematically 
examined  with  regard  to  the  convergence  of  the  time-dependent  numeri¬ 
cal  solutions,  and  the  point  is  made  that  the  initial  conditions  can 
not  be  completely  arbitrary.  Finally,  in  order  tc  learn  more  about 
the  performance  of  the  time-depcr-dent  method,  a  numerical  experiment 
is  conducted  to  examine  the  unsteady  propagation  and  region  of 
influence  of  a  slight  pressure  disturbance  introduced  at  a  point  on 
the  surface  of  a  blunt  body. 
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